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We have recently introduced a discrete model of Lorentzian quantum gravity, 
given as a regularized non-perturbative state sum over simplicial Lorentzian space- 
times, each possessing a unique Wick rotation to Euclidean signature. We investigate 
here the phase structure of the Wick-rotated path integral in three dimensions with 
the aid of computer simulations. After fine-tuning the cosmological constant to 
its critical value, we find a whole range of the gravitational coupling constant ko 
for which the functional integral is dominated by non-degenerate three-dimensional 
space-times. We therefore have a situation in which a well-defined ground state 
of extended geometry is generated dynamically from a non-perturbative state sum 
of fluctuating geometries. Remarkably, its macroscopic scaling properties resemble 
those of a semi-classical spherical universe. Measurements so far indicate that fco 
defines an overall scale in this extended phase, without affecting the physics of the 
continuum limit. These findings provide further evidence that discrete Lorentzian 
gravity is a promising candidate for a non-trivial theory of quantum gravity. 
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1 Introduction 



The aim of discrete approaches to quantum gravity is the non-perturbative construc- 
tion of a quantum theory of gravity as the continuum hmit of a discretized state 
sum, analogous to Feynman's construction of the quantum-mechanical propagator 
of a particle [|1] . Among quantizers of gravity, path- integral formulations have fallen 
into disrepute, both because of the non-renormalizability of the perturbation series 
and because of the unboundedness of the action, which seems to render Euclidean 
approaches (at least formal continuum path integrals and covariantly formulated 
cosmological models) ill-defined. Although this does not necessarily constitute an 
obstacle to the existence of a non-perturbative path integral, it raises the question 
of how such a quantity is to be constructed. 

There is little hope of evaluating the continuum path integral directly, because of 
the complicated functional form of the gravitational action and because of quantum- 
field theoretic divergences, with the ensuing need to regularize in a way compatible 
with the symmetries of the theory. This leads to additional complications, since the 
invariance group of general relativity is the group of space-time diffeomorphisms, 
and not the Poincare group of quantum field theory on a fixed, flat background, so 
that standard regularization methods cannot be applied. 

There are alternative regularization schemes, much in the spirit of Feynman's 
treatment of the non-relativistic particle, formulated on discretized versions of the 
space of all space-time geometries (that is, of space-time metrics modulo diffeo- 
morphisms). One popular class of models is based on applying Regge's idea of 
approximating smooth space-times (M, g^^) by piecewise linear, simplicial manifolds 
in a quantized context. Unfortunately such models have to date produced little con- 
vincing evidence of interesting continuum physics, for reasons that are ultimately 
not well understood (see for a recent review of discrete models in 4d). 

A serious criticism one can level at these models is that they are all formulated for 
positive-definite Euclidean (Riemannian) metrics g^^, and not for physical metrics 
g^°^ of indefinite, Lorentzian signature. This is done for technical reasons, since in 
a concrete regularized formulation one must make sure that the state sum/integral 
Z converges. However, it is important to realize that in a non-perturbative context 

Lor(M) Eu(M) 
Diff(M) Diff(M) 

is an ad hoc substitution: away from a handful of metrics with special symmetries 
(for example, flat Minkowski space), there is no straightforward "Wick rotation" 
9]l°' ^ (or, equivalently, t -it). 

It may well be that the absence of an interesting continuum limit in these sta- 
tistical models of dynamical geometries is related to the absence of any Lorentzian 
structure in their partition functions. This observation has motivated us to construct 
a well-defined discrete quantum gravity model in terms of Lorentzian geometries. A 
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suitable starting point is the method of dynamical triangulations, a variant of the 
quantum Regge calculus program, which has the advantage of being amenable to 
both numerical simulations and analytic treatments. Following Regge's original 
concept of describing "geometry without coordinates" , the "sum over all paths" is 
performed directly over physically inequivalent geometries. Unlike in continuum 
path-integral approaches, there is no need to introduce coordinates and to subse- 
quently gauge-fix them. (In this sense, diffeomorphism-invariance is manifest.) In 
our model, instead of using equilateral Euclidean triangulations, we take the state 
sum over a certain class of Lorentzian Regge manifolds, obtained by gluing together 
a number of simple simplicial Lorentzian building blocks. 

The model has been constructed explicitly in 2, 3 and 4 space-time dimensions 
1^ ^ H, and been solved exactly in d=2 0. Each Lorentzian geometry (or "his- 
tory") has a foliated structure, with a (discrete) proper time t labelling successive 
spatial slices. (Note that this proper time is simply one of the parameters char- 
acterizing the discrete geometries, and not a "gauge choice", since the formalism 
is completely coordinate- invariant from the outset.) In addition, each history has 
a causal structure, induced from the piecewise linear Lorentzian metric structure. 
Each spatial slice is a (d— l)-dimensional equilateral triangulation of Euclidean sim- 
phces, with squared edge lengths ^gp^ce = ct^, and spatial topology changes are not 
allowed (in line with the continuum notion of causality). 

A unique Wick rotation is defined on every Lorentzian history. It maps a given 
triangulation with certain assignments of edge lengths into the same topological 
triangulation, but with the (squared) edge lengths of its time-like links (which in- 
terpolate between the spatial slices) redefined from l^^^^ = —aa^ to /time = 
where a > 0. This leads to an analytic continuation from Lorentzian to Euclidean 
signature of the Regge action in the complex a-plane and gives a precise meaning 
to the map (Q). 

For finite (discrete) volume, the Lorentzian gravity models thus obtained are well- 
defined in the sense of being statistical systems whose transfer matrix is bounded 
and positive. One is then interested in whether they exhibit any critical behaviour 
as functions of the bare coupling constants, leading to continuum theories of quan- 
tum gravity. Secondly, one wants to compare their properties with those of the 
corresponding Euclidean dynamically triangulated models. It should perhaps be 
emphasized at this point that our non-perturbative path-integral method is in prin- 
ciple exact (and not formal). It does involve a specific choice of a "measure", but 
one would expect from universality arguments that the fine details of any choices 
made at the discretized level will not alter the continuum theory. 

As already mentioned above, in d=2 the Lorentzian model is exactly soluble and 
lies in a different universality class from its Euclidean counterpart (which can also be 
solved exactly and is better known as 2d Liouville quantum gravity). Its geometric 
properties are different, which also affects its behaviour when matter is coupled to 
the gravitational degrees of freedom p|, |]. Quite remarkably, and unlike Liouville 
gravity, the coupled system remains well-behaved beyond the so-called "c=l barrier" 
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(in our case, this is equivalent to the number of coupled Ising models exceeding 2). 

These are very interesting results from the point of view of systems of two- 
dimensional random geometries, but our ultimate interest lies in the physical, four- 
dimensional case, and the physics of "general relativity" in 2, 3 and 4 dimensions is 
very different. Where dynamically triangulated Euclidean models seem to go wrong 
in (i > 2 is in the dominance of highly degenerate geometries over their statistical 
ensembles. It is encouraging that one can show the absence of the same type of 
geometries from the Lorentzian ensemble ^, but one could still be worried about 
the occurrence of (less extreme) pathologies. The only way to determine whether 
Lorentzian gravity does indeed solve the problems of the Euclidean approach, is to 
investigate its phase structure in the continuum limit, either by numerical simula- 
tions or by solving it explicitly. 

Before embarking on the physically relevant case of d = A, we will in this paper 
investigate Lorentzian quantum gravity in three dimensions. (Some of the results 



presented here have been announced recently in |1T0| , |TT|.) Apart from being a new 
statistical model of three-dimensional fluctuating geometries, this also has some 
interest from the point of view of quantum gravity proper. Although largely an 
unphysical theory, 3d quantum gravity is an extensively studied system |T^, |T^ . It 
is often invoked as a model system for the full theory, since its classical equations 
resemble in many ways those of general relativity. The big difference from d=4 is the 
fact that there are no propagating physical fleld degrees of freedom. After solving the 
constraints of the theory, only a finite-dimensional phase space remains. Although 
one has not yet been able to make full use of this observation in a configuration 
space path-integral formulation, it suggests that one may still be able to solve 3d 
gravitational models analytically. 

The derivation of the partition function and the explicit construction of 3d 
Lorentzian simplicial space-times was given in In order to make this article 
self-contained, we will summarize the main results below and at the beginning of 
Sec. 3. The Einstein action of a given (smooth) Lorentzian geometry in three di- 
mensions is 

SM = J d^x ^/-detg {R - 2A), (2) 

where G and A denote the gravitational and cosmological constants. For the con- 
tinuous, piece-wise linear geometries employed in our simplicial discretization, we 
use the Regge form of the action p|, |14[, expressed purely in terms of geometric 
(coordinate-independent) data, namely, the geodesic edge lengths of the simplicial 
complexes. 

In order to make the state sum well defined, we analytically continue the Regge 
action associated with each 3d Lorentzian triangulation to Euclidean form by chang- 
ing the length assignments of all time-like edges from l^^^^ = —a^ to = +a^. The 
main aim of this article is to analyze the phase structure of the model defined by 



4 



the path integral obtained after this "Wick-rotation", 



cm ^ cm ' 

where C (T) is the order of the symmetry group of the triangulation T . The "measure 
factor" 1/C(T) appears naturally in the counting of unlabelled triangulations |T5[ . 



As usual in the theory of critical phenomena, we do not expect the detailed choice 
of the measure to affect the continuum limit of the theory, a behaviour that has 
already been corroborated by the 2d models of Euclidean and Lorentzian quantum 
gravity. In (^, the Euclidean three-dimensional Regge action S'e is expressed in 
terms of the total numbers of vertices and tetrahedra, A^o and A^3, according to 

^e(A^o, A^3, T) = -koNo + hN^, (4) 
with the associated dimensionless bare coupling constants 

ko = ^, ks= ""'^ +4^(3k-1) (5) 

(see appendix 1 for a derivation). This form of the action is familiar from former 



work in Euclidean dynamical triangulations [It]. The parameter kg is proportional 
to the bare inverse gravitational coupling constant, while k^ is a combination of 
the bare gravitational and cosmological constants (often referred to - somewhat 
imprecisely - simply as the (bare) cosmological constant). The geodesic "lattice 
spacing" (edge length) is given by a > and kt: = arccos(l/3) is the dihedral angle 
of an equilateral tetrahedron. 

To keep things simple, we are assuming that the spatial slices have the topology 
of two-spheres. In addition, for the convenience of the numerical simulations, we 
are using periodic boundary conditions in the (Euclidean) time direction, unless 
specified otherwise. The sum in (|^) is taken over the set of all causal triangulations 
Tt{S^ X S^) compatible with this topology, and constructed according to the rules 
described in Sec. 2 below. The integer parameter T denotes the total extent in 
(discrete) proper time, i.e. the number of spatial slices of constant t G [0,T].[| 

In this article we explore the phase diagram of three-dimensional discrete Lo- 
rentzian gravity. We are particularly interested in identifying those regions of 
coupling- const ant space where a continuum limit may exist. This is done with the 
help of Monte Carlo simulations of the statistical ensemble defined by (^), combined 
with qualitative analytical arguments. 

The rest of this article is organized as follows: the next section contains some 
general considerations on the behaviour of discrete quantum gravity models un- 
der renormalization. In Sec. 3 we describe the implementation of the Monte Carlo 



^We have slightly changed our notation with respect to §, where the total proper time was 
called t. 
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algorithm on the ensemble Tt{S^ x S'^) of causal 3d geometries, as well as a charac- 
terization of the triangulations and the Monte Carlo moves in terms of dual graphs. 
Our numerical results are presented in Sec. 4. We characterize the different phases 
by measuring various geometric observables, and give a detailed description of the 
geometric properties of the physically interesting "extended" phase. Finally, Sec. 
5 contains a summary and discussion of our results. There are three appendices. 
In App. 1, various identities and derivations for 3d simplicial geometries are col- 
lected, App. 2 contains more details on dual graphs, and App. 3 some technical 
specifications of the Monte Carlo simulation. 

2 Renormalization in discrete quantum gravity 

As we know from the study of Euclidean simplicial quantum gravity |T^, there is 
a well-defined strategy to search for possible continuum limits for the type of dis- 
cretized gravity model we are considering. For each value ko of the bare inverse 
gravitational coupling there is a critical value ^3(^9) of the bare cosmological con- 
stant such that the model is well defined for k^, > k^^ko) and diverges for k^ < k^i^ko). 
One can hope to obtain a continuum limit for k^ ^3(^0) because in this limit the 
expectation value (A*"^) may diverge for suitable powers n. 

This program has been carried out successfully in both Euclidean and Lorentzian 
quantum gravity in d=2 (see [0 for a recent review). It offers a non-perturbative 
field-theoretical definition of 2d quantum gravity where the bare cosmological con- 
stant k2 is additively renormalized according to 

k2 = k^2 + Aa^ (6) 

with A denoting the 2d continuum cosmological constant, and where the critical fcg 
comes entirely from the entropy of the two-dimensional triangulations. We expect 
an analogous additive renormalization of the bare cosmological coupling constant 
^3 in 3d quantum gravity, but in this case k^ depends not only on the entropy of 
the three-dimensional triangulations but also on the gravitational coupling ko, since 
the Einstein action gives a non-trivial weight to each triangulation (contrary to two 
dimensions, where the curvature term is topological). 

Taking the infinite-volume limit of a regularized quantum field theory does not 
necessarily lead to a continuum quantum field theory. For example, the Ising model 
on a infinite two-dimensional lattice will only represent a c = 1/2 conformal field the- 
ory if at the same time the temperature (which plays the role of a coupling constant 
in the theory) is fine-tuned to the critical temperature of the Ising model. Only 
when approaching the critical temperature will the long-range spin fluctuations be- 
come important and allow us to forget about the details of the lattice regularization, 
thereby making contact with continuum physics. 

By contrast, 2d quantum gravity is an example where the infinite- volume limit 
of the regularized theory automatically leads to the continuum theory. This was to 
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be expected as the cosmological coupling is the only coupling constant of the theory, 
and at the same time conjugate to the space-time volume. 

It is not immediately clear what to expect in three-dimensional quantum gravity. 
The classical theory (after gauge-fixing) has no propagating field degrees of freedom, 
but is described by a finite number of (Teichmiiller) parameters, whose number de- 
pends on the topology of the spatial slices. Since in the theory of critical phenomena 
a divergent correlation length and the associated fine-tuning of a coupling constant 
are usually associated with a field degree of freedom, it is tempting to conjecture that 
the situation will be as in two dimensions, namely, that the infinite-volume limit of 
3d Lorentzian quantum gravity (obtained by fine-tuning the cosmological constant) 
coincides with the continuum limit, without the need for further fine-tuning. 

In three dimensions, we must in addition understand which role the gravita- 
tional coupling constant plays in our formulation. In the exponentiated action, it 
multiplies the curvature term J d^x ^ydet g{x) R{x) of the classical Einstein action, 
from which the classical dynamics is derived. Moreover, it is exactly this term that 
gives rise to the non-renormalizability of three-dimensional quantum gravity, when 
one considers perturbation theory around a classical solution. This means that it 
does not make much sense to expand around a given flat background in a conven- 
tional way. Although the underlying quantum theory may not have any divergences 
(since there may not be any propagating field degrees of freedom), we are likely to 
end up with a infinite set of divergent countertermsQ, unless there is some as yet 
undiscovered cancellation mechanism. 

There are well-known ways to circumvent this deadlock in the case of three- 
dimensional gravity, for example, by quantizing in the reduced, finite-dimensional 
phase space, either in a geometric formulation using g^j^y or using gauge-theoretic 
(Chern-Simons) variables [|12|, (However, it should be remembered that even 
classically, the relation between these two "time-full" and "time-less" formulations is 
only partially understood |1^.) How this is reflected in a path-integral quantization 
in terms of geometries is much less clear (we mean here a configuration space path 
integral, rather than a phase space path integral). As far as we know, there is 
not even a generally accepted answer to whether or not the gravitational coupling 
constant should be renormalized. 

If our discretized non-perturbative model possesses a continuum limit, it should 
provide an answer to this question. For example, the presence of a non-trivial 
second-order phase transition at a specific value of would strongly suggest to 
take the continuum limit by fine-tuning to this point, defining in the process the 
renormalization of the gravitational constant. The issue of non-renormalizability 



^One could try to view such a theory as an effective low-energy theory with limited predictive 
power, much in the same way as the non-linear sigma model is used as an effective field theory for 
pion physics, describing some aspects of low-energy QCD. However, this is not what we are after 
in a non-perturbative definition of 3d quantum gravity. We want a theory whose predictions in the 
continuum limit are all independent of the cut-off of the regularized theory, although they may in 
principle depend on a non-perturbatively induced mass scale. 
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could be circumvented if the fixed point was non-trivial, not allowing for a simple 
expansion in geometry. To some extent this is realized in 2 + e-dimensional quantum 



gravity which possesses a non-trivial fixed point [|T9|, |21|. Of course the challenge 
in such a scenario would be to understand how its excitations are related to the 
degrees of freedom, or rather the lack of degrees of freedom, of the classical 3d 
gravity theory. Conversely, if no second- or higher-order transition is present and 
one can still define a continuum limit, it is likely that no renormalization of the 
gravitational constant is necessary. 

In order to find answers to these questions, we will in this article analyze data 
coming from numerical studies of three-dimensional Lorentzian gravity. Attempts 
to solve the model analytically are under way and will be reported elsewhere ||2^ . 
As will be described in the following, our investigation provides evidence that 



a continuum limit exists 



(ii) there is a well-defined ground state which dominates the functional integral 
and thus represents a "background geometry" ; 

(iii) the bare gravitational coupling constant sets a length scale for the geometry, 
but is not renormalized. 



3 Numerical implementation of the model 

Which are the three-dimensional Lorentzian geometries contributing to the state 
sum (|l])? Starting from a sequence of two-dimensional equilateral triangulations, a 
three-dimensional Lorentzian triangulation is obtained by filling the spaces between 
pairs of such spatial slices by three types of tetrahedral building blocks, in such a 
way as to form a simplicial manifold. They are (i) the so-called (3,l)-tetrahedra 
with a triangle in the spatial t-plane and a vertex in the spatial t-|-l-plane; their 
number in any given sandwich [t, t + l] is denoted by A^3i(t), and their total number 
by (ii) the (l,3)-tetrahedra with a vertex in the t-plane and a triangle in the 
t + l-plane; their number in any given sandwich [t,t + 1] is denoted by Nis{t), and 
their total number by and (iii) the (2,2)-tetrahedra with one link in the t-plane 
and another one in the t + l-plane; their number in any given sandwich [t,t + 1] is 
denoted by N22(t), and their total number by N22 (see Fig. p. 

Several of our numerical measurements involve the two-volume of spatial slices. 
In our model there are two natural ways of defining the spatial volume at time t. 
One can define it simply as the number of triangles in the spatial slice of constant 
integer t, 

N^^'\t) = N3,it) = N,,it - I) (7) 
or measure it at half-integer t and define^ 

Nt'^'\t) := N,,{t) + Nuit) + N22{t). (8) 

•^This definition is the simplest one in that it counts the number of building blocks at t — 1/2. 



8 




Figure 1: The three types of tetrahedral building blocks used in discrete 3d 
Lorentzian quantum gravity, and their location with respect to the spatial slices 
of constant integer t. 

In a phase of extended geometry (such that A^i3 ~ N^i ~ N22), both definitions 
should lead to equivalent results. For most purposes, we have found it convenient 
to work with the two- volumes A^2'*'''^''^''(^)- 

We will explore the infinite- volume limit of the ensemble of discrete Lorentzian 
geometries by performing a Monte Carlo simulation where each suggested local 
change of triangulation (a "move") is accepted or rejected according to certain prob- 
abilities depending on the change in the action and the local geometry. 

Our local updating algorithm consists of five basic moves. They change one 
Lorentzian triangulation into another, while preserving the constant-time slice struc- 
ture, as well as the total proper time T. We are confident that this set of moves is 
ergodic in the space of all allowed Lorentzian triangulations at fixed T, although we 
do not as yet have a complete formal proof. Note that all of the moves described 
below will be rejected in the updates if they lead to triangulations where pairs of 
vertices are connected by more than one link or where triplets of vertices belong to 
more than one triangle, since this violates the simplicial manifold property. Let us 
now describe each of the moves in turn: 

(1): Consider two neighbouring triangles in the spatial t-plane. Each of them 
belongs to a tetrahedron above and below that plane. Assume now that both 
the two (3,l)-tetrahedra above and the two (l,3)-tetrahedra below share a 
triangle. Together, the four tetrahedra form a diamond whose tips lie in the 
t—1- and in the t-f-l-plane, and whose intersection with the t-plane is a square. 
The move consists in flipping the link that forms the diagonal of this square 
to the opposite diagonal, accompanied by the corresponding reassignment of 
the tetrahedra constituting the diamond (Fig. 

(2&3): Consider a triangle in the t-plane together with its two neighbouring tetra- 

We could put in individual weights reflecting the actual areas of the cross-sections of the tetrahedra, 
but it would not affect our results below. 
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Figure 2: In a flip move, the four tetrahedra inside a diamond conflguration are 
reassigned by flipping the diagonal of the central square. 



hedra, whose two tips Vt+i and Vt-i lie in the t + 1- and the t — 1-plane. We 
can always insert a vertex Vt at the centre of the triangle and connect it to 
the exterior vertices of this conflguration by adding flve internal links, thus re- 
placing the original two tetrahedra by six (Fig. |^) . The corresponding inverse 
move can only be performed if we can identify a vertex vt of order six (i.e. 
belonging to six tetrahedra), together with two links {vt,Vt-i) and {vt,Vt+i) 
which are both of order three. In this case one can just remove Vt and both 
links {vt,Vt±i), replacing in an obvious way the six tetrahedra by two. 

(4&5): The fourth move can be performed on any conflguration consisting of a pair 
of a (2,2)- and a (3,1)- (or a (1,3)-) tetrahedron having a triangle in common. 
We can remove the triangle (but not its links and vertices) and insert a link 
dual to it, connecting the two vertices which did not belong to the triangle 
(see Fig. H). In this way the original (3,1)- and (2,2)-tetrahedra are replaced 
by one (3,1)- and two (2,2)-tetrahedra, without introducing any changes in 
the two-dimensional spatial slices. The flfth move is the inverse of the fourth 
move, replacing a suitable conflguration of one (3,1)- (or (l,3)-)tetrahedron 
and two adjacent (2,2)-tetrahedra by a pair of a (3,1)- (or (1,3)-) and a (2,2)- 
tetrahedron. 

Note that not all of the local moves preserve the three-volume. We will use a 
standard way of dealing with this situation, developped for dynamically triangulated 
models in dimensions three and four |]16|. This method ensures that the volume of 



the system fluctuates around a prescribed value N^, with a well-deflned range of 
fluctuations. 

In the implementation of the numerical code it is convenient to work not with the 
triangulations and their constituents but with the dual graphs, which are given by 
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Figure 3: Insertion or removal of a vertex in the central triangle transforms two into 
six tetraliedra and vice versa. 




Figure 4: A neighbouring pair of a (3,1)- and a (2,2)-tetrahedron is replaced by 
another (3,l)-tetrahedron (whose tip lies now in the top right-hand corner) and a 
(2,2)-tetrahedron on either of its flanks. 

specific classes of 0^-graphs. Like the triangulations, all of the graphs have a foliated 
structure. This is most naturally associated with half-integer times, because the 
vertices of the dual graph are located at the centres of the tetrahedra of the original 
triangulation. To visualize the geometry of the gluings and the moves in this dual 
language, we adopt a colouring for the dual graphs. A link dual to a triangle 
of a (3,l)-tetrahedron is "blue", and one dual to a triangle of a (l,3)-tetrahedron 
is "red". (We have already anticipated this in Figs. |1]-^ by giving each triangle 
a definite colour; only the links affected by the Monte-Carlo moves are drawn in 
black.) This results in a unique colouring for all links dual to "time- like" triangles 
(lying in between spatial slices), since it is not possible to glue directly a red to a 
blue triangle. This can only be done if the triangles are space-like (i.e. if they are 
both contained in the same slice t = const). The links dual to such triangles are 
therefore double-coloured (Fig. ^ left). We can now construct for each "sandwich", 
i.e. each triangulated space-time slice [t,t+l], a bi-coloured graph (with topology 
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Figure 5: How the tetrahedral building blocks give rise to dual bi-co loured graphs 
(drawn as red and blue horizontal lines) at half-integer t. 



S'^) by projecting all uni-coloured dual links associated with the sandwich to the 
plane at t + 1/2. 

In this way each (3,l)-tetrahedron gives rise to three blue links, sharing a triva- 
lent intersection. Each of the links can end at the centre of either another (3,1)- 
tetrahedron or a (2,2)-tetrahedron, but never at the centre of a (l,3)-tetrahedron. 
(An analogous statement holds for the triplet of red links associated with a (1,3)- 
tetrahedron.) Consequently, each (2,2)-tetrahedron in the sandwich corresponds to 
a four-valent vertex of the dual graph, with alternate colours blue-red-blue-red for 
the incoming links (Fig. |^, right). The end result is a combined red-and-blue graph 
in the t + l/2-plane. Moreover, each such graph occurs in the large-n limit of the 
perturbative expansion of the two-matrix model defined by the partition function 

Z{an, aB, P) = J d(f)R d^^ e"*^ [-i(</'li+0|)+^0|+^</'l+|<A«0B<^H<^s]^ (g) 

where, as usual, the quadratic terms give rise to propagators or links, and the cubic 
and quartic interaction terms correspond exactly to the tri- and four-valent inter- 
sections illustrated in Fig. ||. Note that not all graphs generated by (|]) correspond 
to allowed Lorentzian triangulations, since their duals may violate the 3d simplicial 
manifold constraints. This matrix model (with some additional assumptions about 
universality) can be taken as the starting point for an analytical solution of the 
transfer matrix of simplicial 3d quantum gravity P2 . 



The time-evolution in the dual picture can be thought of as follows. A bi- 
coloured graph at time t + 1/2 consists of two components: a blue 0^-graph dual to 
the triangulation at time t, and a red 0^-graph dual to the triangulation at t+1. The 
way in which the two original spatial triangulations are glued together is encoded 
in the intersection pattern of the "superposition" of the two graphs at time t+1/2. 
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If we view the blue and red trivalent graphs as representing in- and out-states, their 
transition amphtude is a function of the number of topologically inequivalent ways 
of superposing the two graphs (subject to some "dual" manifold constraints - see 
App. 2 for details). 

The five Monte Carlo moves described earlier can also be rephrased in the lan- 
guage of intersecting coloured 0^-graphs, as illustrated in Fig. The diagrams 
appearing in Fig. ^ contain all dual links affected by a given move. Since the moves 
1, 2 and 3 are symmetric with respect to the plane t = const, a change in one of the 
trivalent graphs is always accompanied by an equivalent change of its mirror image 
of the opposite colour. Moves 4 and 5 assume a particularly simple form: one link 
of a given colour is "dragged across" a vertex of the opposite colour. (Note that 
in the graphical representation of this particular move, the cubic vertices represent 
(3,1)- or (l,3)-tetrahedra and the blue-red crossings (2,2)-tetrahedra.) Details about 
the numerical implementation in terms of this dual picture (including lattice sizes, 
update efficiency, number of sweeps etc.) can be found in Apps. 2 and 3. 



4 Numerical results 

Having presented our numerical set-up, we are now in a position to extract a number 
of physical properties of the Lorentzian model. We will first investigate the phase 
diagram of the regularized theory, and try to understand which of the continuum- 
limit scenarios outlined in Sec. 2 is realized. We will then analyze the geometry 
characterizing the different phases. Since we have a distinguished (and coordinate- 
invariant) notion of proper time t, we can extract invariant information of the system 
by studying correlation functions in t. This will be done by measuring distributions 
of spatial slice-volumes N2^^^'^^ as a function of the total proper time T and cor- 
relators (A^2'^'''^'^^''(^i)^2'*'''^'^^''(^2)) between spatial volumes, as well as the intrinsic 
Hausdorff dimension (f^ of a typical spatial slice. 



4.1 The phase diagram 

In order to explore the phase diagram of the regularized Lorentzian model we must 
find an order parameter, and explore how it changes with the coupling constant, 
in this case ko. We have found that the ratio between the total number A'"22 of 
(2,2)-tetrahedra and the total space-time volume N^, 

N22 N22 , ^ 

T= — ~ , 10 

serves as an efficient order parameter. We shall not be concerned with a continuum 
interpretation of this parameter (which is not obvious) since we will go on to show 
that no continuum physics is associated with the transition we observe as a function 
of r. In Figs. |^ and § we show the ratio r as a function of fco for two different 
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Figure 6: Graphical illustration of the Monte Carlo moves in terms of dual bi- 
coloured graphs. In moves 1 and 2&3, dual links from two adjacent sandwiches are 
affected. Moves 4&5 take place within a given [t, t + l]-sandwich. 

types of space-time configurations. In Fig. |^ all geometries have 64 spatial slices 
(T = 64), with total space-time volumes = 16, 000 and A^3 = 64, 000. One observes 
a rapid drop to zero of T{ko) around ko ^ 6.64. Increasing N3, the drop becomes 
a jump, characteristic for a first-order phase transition. A detailed study of the 
neighbourhood of = 6.64 reveals a (weak) hysteresis as one performs a cycle, 
moving above and below the critical value /cg, again as expected in a first-order 
transition. The location of the phase transition depends weakly on the total length 
T in time-direction, and the jump in r becomes more pronounced for larger volumes. 
Fig. ^ is an illustration of the extreme case T = 1, where we have chosen free instead 
of periodic boundary conditions (the triangulations at the initial time t = l and the 
final time t = 2 are allowed to fluctuate freely), so that the topology of space-time is 
changed to x [0,1].^ 

^In this situation the critical point is changed from fco = 6.64 to fco = 6.42. 
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Figure 7: The order parameter r = N22/N3 for configurations with T = 64, and 
A^'3 = 16,000 (long dark curve) and 64,000 (short hght curve), plotted as a function 
of ko- The curve is a linear interpolation between data points. (Error bars smaller 
than width of curve.) 

As can be read off from Fig. 0, for ko > ^ 6.64 only a minimal number of 
(2,2)-tetrahedra is present. This can be understood by rewriting the action (^) to 
make the dependence on N22 explicit. In appendix 1 we derive 



which shows that for fixed and T (and positive gravitational coupling ko) a 
minimal N22 corresponds to a minimum of the Euclidean action. 

The entropy of configurations with N22 different from its minimal value will in 
general ensure that the ratio r is different from zero, even when A^3 — > 00. However, 
since the number of such configurations for fixed A^3 grows at most exponentially 
with A^3, this leaves the possibility that for sufficiently large ko the term e~^''^22/4 
will dominate over the entropy contribution and trigger a phase transition to a phase 
with only a minimal number of (2,2)-tetrahedra, such that r = in the continuum 
limit. 

The physics of this phase can be readily understood. In terms of the matrix 
model, a situation with no (2,2)-tetrahedra corresponds to choosing the coupling 
constant /3 = in (^), thus reducing the model to a product of two independent 
0^-matrix models. Since a 0^-matrix model at its critical point describes two- 
dimensional Euclidean quantum gravity, the matrix model analogy strongly suggests 
that the Lorentzian 3d model for ko > k^ can be viewed as a product of uncoupled 
2d gravity models. This conclusion seems to be corroborated by our numerical re- 
sults. Fig. I is a typical "snapshot" of a space-time geometry, taken during the 
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Figure 8: The order parameter r = N22/N3 for configurations with T = 1 and 
= 16, 000, but with free boundary conditions, plotted as a function of fco- The 
curve is a hnear interpolation between data points. (Error bars smaller than width 
of curve.) 



computer simulations. The spatial volume N2^^^'^\t) is shown as a function of 
the time t. Apparently it can change from essentially zero to a "macroscopic" size 
in a single time-step, which implies that there cannot be any correlations between 
slices separated by a few time-steps. A different measurement of the correlation 
between successive spatial volumes is depicted in Fig. |l^. We have again chosen 
T = 1 and free boundary conditions, in order to have spatial slices of a reasonably 
large size, but a qualitatively similar behaviour is expected for T > 1 too. We show 
the distribution of the (normalized) difference 

..^ ^ |Arj-)(t = 2)-Arj-)(t=l)| ^ |ArW(2) - ivW(l)l .^2) 

of the spatial volumes of the initial and final slice. For /cq less than the critical 
/cq ~ 6.42, f is peaked around zero. The peak becomes flatter as /cq approaches its 
critical value and immediately beyond k^, the distribution approximates a (5-function 
around f = 1 . 

This last result can be understood as follows. We know from the simulations that 
the number of (2,2)-tetrahedra drops to a minimum beyond the critical point k^. To 
first approximation, these tetrahedra therefore do not contribute to the entropy in 
that region. Moreover, since a minimal set of (2,2)-tetrahedra can basically inter- 
polate between any pair of "incoming" (3,1)- and "outgoing" (l,3)-configurations, 
the combinatorics is governed by the separate countings of those configurations, 
subject only to an overall volume constraint N^i + Ni^, = N3 — N22 ^ = const. 
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Figure 9: Monte Carlo snapshot of the distribution of spatial volumes N2^^^'^\t), 
for T = 32, A^3 = 16,000 and A;o = 6.7 (that is, above the critical value fco = 6.64). The 
volumes are plotted symmetrically about a central horizontal t-axis. 



Individually, the configurations at t = 1 and t = 2 are simply 2d Euclidean triangula- 
tions, whose number for a given spatial volume N^^^ is known to be proportional to 
e'^^2''^(Ar^^))-5/2_ pi-oin f (1 ^ 2) = |1 - 2iVi3/iV3|, and taking into account that the 
minimal interpolating Ai"22-configuration can be inserted anywhere in the incoming 
and outgoing configurations, one finds 



#(iV3i ^ iVi3) ~ e^^-Ar37e^(^3-iV30(iV3 _ n^,)-- ~ (1 - f'} 



(13) 



for the combined entropy at fixed volume N^. In agreement with Fig. |T0[ it shows 
that the most likely configurations are those where the entire 3d volume is con- 
centrated at one of the slices, that is, either N31 ~ 0, A^i3 ~ A^; 



3 or vice versa. 



The situation is very different in the phase with < /cq. Our measurements 
of f at T = 1 are an indication that also in general in this phase the volumes of 
successive spatial slices will be strongly coupled, i.e. their volume difference will be 
small. Since this behaviour is not favoured by the action, the prevalence of such 
configurations must have to do with the presence of the (2,2)-tetrahedra and their 
associated entropy (i.e. the number of ways they can combine with each other and 
with other tetrahedra to form interpolating "sandwiches"). This observation turns 
out to be of great importance, since it seems to lie at the heart of the emergence of 
extended geometries in this phase, which will be described in more detail in the next 
section. In broad terms, the phase is characterized by < r < 1. In principle there 
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Figure 10: The probability distribution P(f) of f = |A^^'^(2) - nP{1)\/N3 for ko = 
6.0, 6.1, 6.2, 6.3, 6.4, 6.42 (highest to lowest peak around f = 0; the value = 6.42 
corresponds to the phase transition point) and 6.44 (distribution peaked around 1); 
for total volume A'"3 = 8,000 and free boundary conditions. 



may be another phase transition at some smaller (possibly negative) value ko, such 
that r = l for /co < ko (with no (3,1)- or (l,3)-tetrahedra surviving in the continuum 
limit). Indeed, for fixed and negative ko the action ([TTl) has a minimum for 
r ~ 1. (Configurations with r — > 1 for A^3 — ^ cx) can actually be realized.) 

Whether or not the system will undergo a phase transition for sufficiently small 
ko will depend on the balance between action and entropy, which cannot be deter- 
mined by simple qualitative considerations. We have not investigated this region of 
the coupling constant space further, given the limited importance of negative grav- 
itational coupling constants from a quantum gravity viewpoint, and the fact that 
our computer algorithm is not efficient at small ko. 

In summary, we have arrived at the following tentative description of the phase 
diagram of 3d Lorentzian quantum gravity: the bare inverse gravitational coupling 
constant has two critical values, fcg fcg (possibly with fcg = — oo). For ko > 
/cq the model describes the fluctuations of an uncorrelated set of two-dimensional 
spatial geometries and has nothing to do with a three-dimensional theory of gravity. 
Also for ko < fco the space-time geometry degenerates, since the spatial slices at 
integer t completely disappear from the theory. These two "extreme" regions of 
the phase diagram can be regarded as artifacts of our particular way of setting up 
the discretized theory. They may be seen as remnants of the phases of degenerate 



geometries observed previously in 3d Euchdean quantum gravity |16|. However 
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unlike the Euclidean theory, Lorentzian gravity possesses a large region fcg < /cq < 
/cq of coupling constant space where the quantum geometry is extended and well- 
behaved, and whose description we shall turn to next. 



4.2 The phase of extended geometry 

Let us now analyze the structure of the phase of intermediate gravitational coupling, 
fco < ^0 < ^0; where all types of tetrahedral building blocks contribute non-trivially. 
Quite remarkably, and unlike in the phase where ko > we observe here the emer- 
gence of well-defined three-dimensional configurations. Fig. |TT| shows a snapshot of 
a typical geometry at ko = 5.0, consisting of 16,000 tetrahedra, for T = 32. (As in 
the previous Fig. we plot - symmetrically around an arbitrary axis - the spatial 
volume N^'^^^^\t) function of t.) Following the computer-time history of this 
extended object, it is clear that although it does indeed fluctuate, the fluctuations 
take place around a three-dimensional object of well-defined linear extension.^] The 
emergence of a ground state of extended quantum geometry is a highly non-trivial 
property of the Lorentzian model, since we have at no stage put in a preferred back- 
ground geometry by hand. No structures of this kind have ever been observed in 
dynamically triangulated models of Euclidean quantum gravity. It underscores the 
fact that the Lorentzian models are genuinely different and affirms our conjecture 
that in c? > 3 they are less pathological than their Euclidean counterparts. 



For a fixed kg and N^, an overall "spherical" shape as shown in Fig. 11 is found 



only for sufficiently large T. At small T, one observes a uniform distribution of 
spatial volumes iVg (t) as a function of t. As T increases, the bulk of the space- 
time volume "condenses" into a region with a well-defined extent T„ < T in time- 
direction (which we will call the universe), leaving over a thin stalk of minimal spatial 
radius everywhere else along the t-axis. We will from now on choose T sufficiently 
large, so that T > for all volumes under consideration. We are interested in the 
"cosmological" properties of this extended universe, i.e. its geometric properties at 
large scales. Our data for the scaling of the time extent T„ of the universe as a 
function of the total volume are consistent with 

Tu ~ Nl'^. (14) 

(We discuss below how a quantitative measure of is obtained.) Similarly, by an 
independent measurement of the volumes N^f^^^"^^ (t) of spatial slices that lie within 
the universe, we have found a scaling behaviour consistent with 

The relations ([l^ and (|15D support an interpretation of the universe as a genuinely 
three-dimensional object. There is of course no a priori reason that a ground state 

^ A trivial mode of fluctuations are the translations in time-direction. Due to the periodicity of 
the boundary condition the "centre of volume" of the extended configuration performs a random 
walk in the t-direction. 
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Figure 11: Monte Carlo snapshot of the distribution N2^^^'^\t) of spatial volumes, 
for T = 32, Ai'3 = 16,000 and /co = 5.0 (that is, below the critical value fcQ = 6.64). The 
volumes are plotted symmetrically about a central horizontal t-axis. 



in a non-perturbative theory of quantum gravity (if it exists) should bear any resem- 
blance with a (semi-) classical geometry. Let us explain briefly how such geometries 
(and, more specifically, classical solutions) might still make an appearance in this 
context. 

For the simplest compact space-time topology, the solution to the classical Ein- 
stein equations with Euclidean signature and a positive cosmological constant A is 
the round three-sphere (of constant positive scalar curvature) with radius Rs^ ~ 
A~^/^. Solving the equations with the constraint of fixed space-time volume V is 
equivalent to introducing an effective cosmological constant Aefj ~ V~'^^^. We are 
not aware of a classical solution with positive A and topology x S*^ (which is the 
topology used in our simulations). However, for our purposes we can "adapt" the 
S^-solution to this topology by cutting away two small open balls at two opposite 
points of an ^^-configuration with radius i?53 ~ ^e//^ ~ V^^^ and attaching a thin 
cylinder / x S*^ (with spatial radius at the cut-off scale) to the holes. This will pro- 
duce a geometry of the kind shown in Fig. |ll|. Although it is not strictly speaking 
a solution to Einstein's equations, it is "almost as good" from the point of view of 
the path integral, since - independent of its metric properties - the contribution of 
the stalk to the action is negligible (because it does not grow proportionally to the 
three- volume) . 

Suppose for the moment that the round S^-solution corresponded to a (local) 
minimum of the action. Then the singular "solution" of topology S*^ x constructed 
above would also represent a (local) minimum of the action, and would therefore be 
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Figure 12: The correlator C(A) with T = QA and A^3 = 32, 000. Dots are the measured 
values (error bars less than dots), and the curve is fitted from the sphere solution 
described in the text. 



as relevant as the 5'^-solution in the quantum theory. Unfortunately, the argument 
is not quite as simple, because the classical continuum Einstein action is unbounded 
from below, due to the presence of a kinetic term of the "wrong" sign, coming from 
the conformal mode of the metric. However, since the conformal mode is not a 
propagating degree of freedom in either classical general relativity or in canonical 
quantizations, it should not cause any problems in a correct, non-perturbative path- 
integral quantization of gravity, not even in the Euclidean sectoifj. In such a quantum 
theory the effective action should be bounded from below and semi-classical saddle- 
point considerations of the kind made above may again be appropriate. 
We have measured the correlation function 

ci^) = ^,f2(^t'^'^mt'^'\t+A)) (16) 
t=i 

as a function of the displacement A to determine the scaling of Tu with the space- 
time volume N^. This correlator has the advantage of being translation-invariant 
in t and allows for a precise measurement by averaging over many independent 
configurations. From the typical shape of the space-time configurations we expect 



C(A) to be of the order of the spatial cut-off if A > 2Tu. Fig. |12| illustrates the result 
of our measurements of C(A), with the dots representing the measured values. The 
theoretical curve to which we are fitting corresponds to the "fake sphere" described 

^The continuum path integral in proper-time gauge is discussed in p3|. 
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above, with the radius of the 5"^ and the spatial radius of the thin cyhnder attached 
to it as free parameters. For this "spherical" geometry we then perform the integral 
(the sum) in (|16]), without the average (■). As is evident from Fig. 0, the volume 



distribution associated with this fixed geometry gives a rather good fit to our data. 
This provides some evidence that we can ignore the quantum average implied by (■), 
and that our universes behave semi-classically, at least as far as their macroscopic 
geometric properties are concerned. We should mention that our "S''^-solution" is 
not singled out uniquely, since the choice of a Gaussian shape in the t-direction gives 
a fit of comparable quality. 

For various space-time volumes (typically 8, 16, 32 and 64k) we have deter- 
mined the radius Rsa of from the fits to the measured C(A). From this, we have 
finally found a = 0.34 ± 0.02 as the best exponent in the scaling relation 

Rs.m = N^. (17) 

The same value is obtained using other ways to extract T^, lending additional sup- 
port to the three-dimensional nature of our universes. 

We will now take a closer look at the geometry of the two-dimensional spatial 
slices. If they could be described as typical triangulated surfaces in 2d Euclidean 
quantum gravity, they would not behave like smooth 2d geometries (when described 
in terms of geodesic distances), but rather like fractal spaces with Hausdorff dimen- 



sion = 4 |2J, By contrast, typical space-time surfaces contributing to the 
path integral of 2d Lorentzian quantum gravity can be viewed as two-dimensional, 
as shown in 0, ||, |^. 

The spatial slices at constant integer t are obviously Euclidean in nature, but it 
is not immediately clear how they will behave, since they appear as part of a larger 
foliated space-time geometry, and are coupled to each other in a non-trivial way. We 
have tried to extract the Hausdorff dimension d^^ of the spatial slices lying inside 
the spherical universe, using the geodesic distance inherited from the 3d geometry, 
and employing techniques developed in the context of 2d dynamically triangulated 
Euclidean quantum gravity ^ . Unfortunately, the quality of our measurements 
is not very satisfactory, since the spatial volumes A^2*'* (t) of the individual slices are 
rather small (typically of the order of Ik). 

One can obtain better data by using simulations with small T (so that T < T„ 
and no universe can form), but it is not entirely clear whether this will leave the 
spatial Hausdorff dimension unchanged. Our measurements for small T point to a 
value around d^^ = 3.4 ± 0.4 (the measurements for larger T are compatible with 
this value, but their error bars are considerably larger). If our experience with the 
2d Euclidean gravity simulations is anything to go by, this probably implies (i^ = 4, 
but so far this has to remain merely a conjecture. At any rate, these somewhat 
preliminary results highlight that the detailed, microscopic geometry of the universe 
may be rather complicated, although its macroscopic properties resemble that of a 
semi-classical object. 
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Attempts to measure the Hausdorff dimension du oi the entire space-time (as 
opposed to that of individual spatial slices) have not yet led to unambiguous re- 
sults. One wants to confine the measurement to the spherical universe, where again 
one runs into difficulties because of its relatively small radius. In addition, one 
needs a dynamical definition of where the universe begins and ends (along the t- 
direction), and must make sure that the result is independent of the particular 
prescription adopted. From the limited data collected (using the geodesic link or 
dual link distance, in the sense in which this notion is usually defined in dynamical 
triangulations) we conclude that the Hausdorff dimension is most likely larger than 
three. 

Another important result concerns the relation between the geometries of differ- 
ent fco, in the phase where fco < ^q. In the numerical simulations we have observed 
the following: 

(i) the distributions as functions of t can be made to coincide for different fco by 
rescaling the time, t — >• fti{ko)t or alternatively at — > fti{ko)at, where at is 
the link length in time direction. This is illustrated by the N2^^^^'^^ -N2^^^'^^ 



correlator C(A), Fig. |T3|, where we show both the actual and the rescaled 
distributions. 

(ii) The distributions measured in the spatial slices from inside the universe can 
be made to coincide for different k^ by rescaling the spatial link distance a^ — > 
fsp{ko)as, where is the length of the spatial links. This is illustrated in 
Fig. |l^ for the distributions of 2d volumes S{1) of spatial spherical shells of 
(link) radius /, measured for various values of ko. (The shell volume S{1) is 
obtained by counting the number of vertices separated from a given vertex v 
by a minimal link distance /. Note that this spherical shell is precisely what 
is measured to determine the Hausdorff dimension c?^ of the spatial slices.) 

(iii) Within the numerical accuracy we find that /ti(/co) = /sp(^o)- In fact, the 
rescaling of the correlator C{A) (Fig. was obtained by simply using the 



values fsp{ko) (see Table 1) determined from the fit S{1) (rather than by finding 
the best value for fti{ko)). 





/sp(^o) 


3.0 


0.84 ± .02 


4.0 


0.91 ± .02 


5.0 


1.00 


6.0 


1.23 ± .03 



Table 1: The spatial scaling factor fsp{ko), extracted from the distributions S{1) for 
various values of ko. 
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Figure 13: The correlator C(A), eq. ( |I6D with T = 32 and Ns = 16,000 measured 
for ko = 3.0, 4.0, 5.0 and 6.0 (lowest to highest peaks). The upper figure shows 
the actual distribution, the lower one the rescaled distributions for positive A (it is 
symmetric in A). The variable x is a rescaled version of the time distance A. 

On the basis of these correlator measurements we conjecture that the value of 
the bare inverse gravitational coupling ko Gjfcg; i^Grely sets the overall length scale 
of the universe, and otherwise does not affect the physics of the model. 

We should point out that the average total integrated curvature is not indepen- 
dent of ko. Subtracting the cosmological term from the action ([3l|), one finds 




where the parameter r (defined in eq. (p!o|) ) now depends dynamically on ko through 
the ensemble average. Comparing with our measured curve for r in Fig. ^ one can 
read off that the total curvature vanishes around ko ^ 5.0. For smaller ko, it becomes 
negative and for larger ko positive. Nevertheless, in line with our conjecture above 
we expect the curvature-curvature correlators to follow the pattern of the already 
measured correlators (i.e. to observe a simple fco-dependent scaling behaviour), but 
this remains to be verified. 



24 



0.15- 



. 125 - 



0.1- 



. 075 - 



0.05 - 



. 025 - 




5 10 15 20 



Figure 14: The 2d volume S{x) of spatial spherical shells, measured only on slices 
inside the spherical universe, for various values of ko and rescaled. The variable x 
is a rescaled version of the radius I. 



5 Summary and discussion 

In this paper, we have analyzed the phase structure of the discretized model of three- 
dimensional Lorentzian gravity defined in with the help of computer simulations. 



The phase diagram. Fig. should be compared with that of the Euclidean theory, 
depicted in Fig. 16. Although the overall phase structure is similar, with a first-order 
transition at some intermediate value fcg, the quantum geometries of the phases on 
either side of the transition are very different in both cases, as indicated in the 
drawings. 



In the Euclidean case |]I6[, one finds a "crumpled" phase at small ko, dominated 
by configurations of very large Hausdorff dimension d ^ oo (these are simplicial 
manifolds where roughly speaking any two vertices are a minimal distance apart). 
Above the first-order transition at k^, the system is in a branched-polymer phase 
of highly branched geometries (with a fractal dimension dn = 2). Unfortunately, 
neither of these phases seems to have a ground state that resembles an extended 
geometry of dimension d > 3. 



Another approach to Euclidean gravity was advocated in p8[ within the quantum 



Regge calculus program. The phase structure found in the numerical simulations 



of this model resembles those of Figs. 15 and at least superficially in exhibit 



ing a "rough phase" for small and a "smooth phase" for large Newton's constant. 
However, this model is inequivalent to the dynamically triangulated models we have 
been discussing, since in the Regge approach these two phases seem to be separated 
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Figure 15: The phase diagram of Lorentzian dynamical triangulations in three di- 
mensions. 

by a second-order phase transition and associated divergent curvature fluctuations, 
indicating the presence of propagating field degrees of freedom (cf. our general dis- 
cussion in Sec. 2). How this can be related to the topological character of 3d quantum 
gravity manifest in canonical treatments of the theory is unclear. 

The situation in Lorentzian dynamically triangulated gravity is completely dif- 
ferent. Although we find a weak "remnant" of the Euclidean degeneracy for ko > k^, 
where space-time decouples into a sequence of uncorrelated two-dimensional slices, 
there is a whole phase below k^ where the geometry is extended, with macroscopic 
scaling properties characteric of a three-dimensional universe. Quite remarkably, 
this is an example in three dimensions of the emergence of a well-defined ground 
state of geometry in a non-perturbative state sum for gravity. In a continuum lan- 
guage, this is the ground state of an effective action, where entropy contributions 
(in other words, the measure) play a crucial role. Apparently in our model these 
contributions are such that they outbalance potential conformal divergences coming 
from the Wick- rotated action (otherwise a well-defined ground state could not exist). 
From the evidence gathered so far, the physics in this extended phase is independent 
of the precise value of the bare gravitational coupling ko. In the correlation functions 
we have measured, ko merely serves to set an effective overall length scale. 

As argued in Sec. 2, these findings strongly favour a situation where the gravita- 
tional coupling is not renormalized, and no fine-tuning of ko is needed to approach 
the continuum limit. This limit coincides automatically with the infinite (lattice) 
volume limit, which we obtained by fine-tuning the cosmological coupling constant. 
Continuum physics is then extracted by taking the limit as A^3 — > oo and a — > 0, 
while keeping the three-volume Kont '■— N^a^ constant. In this setting, no genuine 
field degree of freedom is present since there is no divergent correlation length asso- 
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Figure 16: The phase diagram of Euchdean dynamical triangulations in three di- 



mensions. 



ciated with fine-tuning to the critical point of a second-order phase transition. 

As a result of our investigations, we have good reasons to believe that 3d Lo- 
rentzian quantum gravity, as defined through our discrete, dynamically triangulated 
model, exists as a continuum theory. Since so far this theory is not given in an 
explicit analytical form, the question arises of how to make contact with already 
existing quantizations of three-dimensional gravity. 

At least in spirit, our formulation is related to canonical approaches using the 
trace of the extrinsic curvature as a time variable, the so-called "York time" , with 
a conjugate Hamiltonian determining the time evolution of the system [p9| , |1^ . 
However, this approach only works for genus g > 1, and the only case where the 
canonical quantum theory and the Hamiltonian operator are reasonably explicit is 
g = l, where the spatial slices have torus topology. We are not aware of any quantum 
observables that have been calculated in the case of spherical slices, and which we 
could try to compare to. (For S'^-slices, there are no non-contractible holonomies 
and the reduced phase space is zero-dimensional.) We could in principle repeat the 
simulations for toroidal spatial slices, although the finite-size effects will be larger 
for this more complicated topology (and for T = 32, 64, our spatial slices are rather 
small) . 

Similarly, it is in principle straightforward to enlarge the Lorentzian model to 
include coupling to matter fields. This has already been done in two-dimensional 
Lorentzian quantum gravity, with a clear motivation in mind, namely, to understand 
the status of the c = 1 barrier in general 2d gravity models. We showed that this 
disease of 2d Euclidean quantum gravity can be avoided by working with Lorentzian 
geometries 0. 

It would provide a strong incentive for considering either of these generalizations 
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if there were definite predictions from continuum formulations of 3d quantum gravity 
with and without matter for observables measurable in the computer simulations 
(for example, correlation functions of the type we have been studying). 

A next important step in our analysis of 3d Lorentzian quantum gravity will be 
the derivation of the explicit form of the Hamiltonian in the continuum limit. We 
can in principle obtain the matrix elements of the transfer matrix T between two 
successive triangulated two-geometries gi, as the solution to a combinatorial problem 
22| . Expanding the matrix elements according to 



{g2\T\9i) = {g2\ e-'^^ \g,) = {g^l (1 - aH + 0{a')) \g,), (19) 

one can extract the Hamiltonian operator H. A similar approach was successful in 
2d Lorentzian gravity, where the regularized transfer matrix could be calculated, 
and its continuum limit taken in a straightforward wayQ. The resulting Hamiltonian 
agreed with the one obtained by continuum formal manipulations in proper-time 
gauge PU|, showing that the educated guesses made in this paper were justified. 



This calculation can be generalized to our 3d Lorentzian gravity model, but the 
matrix-model methods will probably only work in the case of a spherical spatial 
topology. A direct comparison with canonical quantum gravity would then require 
a canonical continuum quantization in proper-time gauge, with spatial topology S"^. 

Let us conclude by pointing out an interesting consequence of our arguments 
that would follow if the second critical point (whose value we did not measure) 
was negative. In this case, the theory with bare coupling constant /cq = would lie 
inside the extended phase. This implies that we could start from a discretized gravity 
action with the cosmological term alone, and still obtain the same continuum theory. 
This may seem a radical suggestion, because the classical theory of the action 



S = A d^x ^/det g{x) (20) 



is trivial (it does not contain any time-derivatives). However, there is no logical 
contradiction, since further non-trivial contributions to the (effective) action can be 
generated through the non-perturbative evaluation of the path integral. In fact, 2d 
Euclidean quantum gravity provides a good illustration of this mechanism. There 
the action is given by eq. (I^)' but substituting d?x — > d?x. Nevertheless, the 
effective quantum theory in conformal gauge is described by the highly non-trivial 
quantum Liouville theory. 

What have we learned from our exploration about our ultimate goal, the con- 
struction of quantum gravity in four dimensions? We invented the discrete Lorentz- 
ian models in the hope that they may lead to a better description of physical four- 
dimensional space-time, which after all has a Lorentzian signature. We also conjec- 
tured in that in the continuum limit the causality constraints imposed on each 
geometry in the state sum may lead to a suppression of the degenerate phases of 

^The formula for the Hamiltonian in contains a typo; see [ p^ for the correct expression. 
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highly fractal geometry found in the Euclidean models for d > 3. From the evidence 
presented in this work, this is indeed what happens in three dimensions. Moreover, 
we saw the emergence of a ground state of extended three-dimensional geometry 
in the Lorentzian case. As already observed in d=2, also in three dimensions the 
Euclidean and Lorentzian models correspond to completely different continuum the- 
ories, re-iterating that these two "sectors" of the gravitational quantum theory are 
not related by a simple analytic continuation in time [^, ^ . 

We are very encouraged by these results, since they indicate that also in d=4 
completely different geometries will dominate the Wick-rotated path integral, com- 
pared with the Euclidean theory. The physics that the four- dimensional model 
should describe, if it were to lead to a non-perturbative theory of quantum gravity, 
must of course be very different from that found in two and three dimensions. In 
particular, the critical behaviour of the regularized theory should reflect the presence 
of physical, propagating field degrees of freedom. In the context of the statistical 
models we are considering, the simplest realization would be in terms of a second- 
order phase transition. This possibility is apparently not realized in the dynamically 
triangulated Euclidean gravity models. However, there is by now plenty of evidence 
that the Lorentzian model for quantum gravity defined in ||^, ^ is sufficiently differ- 
ent to make it a new, promising candidate for a non-trivial non-perturbative theory 
of quantum gravity in four dimensions. 
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Appendix 1 

In this appendix we collect some formulas for dynamically triangulated three-geo- 
metries, which were used in deriving various forms of the discrete Einstein action in 
the main text. We will work in the Euclidean sector of the theory, and for simplicity 
choose all tetrahedra to be equilateral (that is, a=—l and /space = /time = a > 0). 
The curvature of a 3d piecewise linear manifold is concentrated at its links. The 
contribution to the total curvature associated with each link / is given by the link 
length a times the deficit angle 

61 = 211 - y^^9i, 

where the sum is taken over all tetrahedra cxj, i = 1, . . . ,o{l), sharing the link /, 
and 6i is the dihedral angle associated with the z'th tetrahedron. For an equilateral 
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three-complex, all dihedral angles are identical, 



9 = arccos = kit, (21) 
and the curvature term of the Einstein action becomes 

- f d^xv^det g{x) R{x) — >"^aSi = 2na{Ni - 3nNs), (22) 

where we have used that for a closed three-dimensional triangulation 

5^1 = iVi, ^o(/) = 6iV3. (23) 
I I 

Taking into account that the 3-volume of an equilateral tetrahedron is given by 
a^/6\/2, we obtain the discretized form of the Euclidean Einstein action ^ 

S, = _^(;V.-3.7V3) + ^;V, (24) 

where in a slight abuse of language we continue to use G and A to denote the bare 
gravitational and cosmological couplings. We can substitute A^i by the number Nq 
of vertices, using the identity Ni = N3 + Nq, which can be derived from the vanishing 
of the Euler number for any closed 3d manifold, 

^ = - iVi + iV2 - iVa = 0, (25) 

together with the relation N2 = 2N^ (any triangle is shared by two tetrahedra and 
any tetrahedron has four triangles). Substituting this into (P^), we obtain the action 
used in Sec. 1, 

Se = -hNo + hNs, (26) 
with the coupling constants given by 

a , a^A a , . , _, 

fco = — , fcg = ^ — + — (3fi;-l). (27) 



In the numerical investigation of Sec. 3 we discussed the dependence of the action 
on the total number N22 of (2,2)-tetrahedra. This can be made explicit by rewriting 
Nq as a function of and N22- For periodic boundary conditions in the t-direction, 
the total numbers of (3,1)- and (l,3)-tetrahedra are the same, and we have 

iVl3=iV31 = ^(iV3-iV22). (28) 
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Next, we need some identities for the spatial slices at constant integer t. Because 
the slices are topologically two-spheres, the number of vertices in a slice is 

iVo(t) = ^iV3i(t) + 2 = ^iVi3(t - 1) + 2. (29) 
Summing this equation over all t and using (|2^) yields 

T 

No = J] Aro(t) = 2T + ^(iVg - (30) 
t=\ 

and therefore 

Se = + ih - j)iV3 - 2koT, (31) 

which is the form of the action used in Sec. 3. 



Appendix 2 

As discussed in Sec. 3, one can describe the 3d Lorentzian geometries in terms of 
dual graphs, naturally associated with each plane of constant half-integer t. They 
decompose into two cubic graphs of different colour (whose trivalent vertices corre- 
spond to the (1,3)- and (3,l)-tetrahedra of the original triangulation, and which may 
be thought of as the in- and out-states of the transfer matrix. Red and blue lines 
cross at four-valent vertices, corresponding to the (2,2)-tetrahedra of the original 
lattice. 

The planarity of this structure (i.e. the fact that the subgraphs have topology 
S"^) is easily implemented in the program by representing the one- dimensional lines 
of the red and blue graphs as double lines with opposite orientation, as one can do 
in the large-n matrix model (^. In this way one obtains closed loops of oriented 
coloured (single) lines which are dual to the vertices at times t and t+1 of the original 
lattice. The 2d spherical surface may thus be thought of as being covered by (either 
red or blue) patches enclosed inside the loops. 

In the numerical simulations we take care that the triangulations are 3d simplicial 
manifolds to start with, and we accept only Monte Carlo moves which preserve this 
property. In terms of the original triangulation, this means that we only allow the 
creation of configurations where any two vertices can be shared by at most one link, 
any three vertices can be shared by at most one triangle, and any four vertices by 
at most one tetrahedron. 

In terms of the dual graphs, this implies two types of restrictions, the first of 
which have a transparent interpretation in the matrix model: they constrain the 
individual trivalent graphs to have neither tadpoles nor self-energy subdiagrams. 
This ensures that they are regular 2d simplicial manifolds with spherical topology. 
The remaining constraints restrict the ways in which the two coloured graphs are 
allowed to intersect each other. Requiring the absence of double links between 
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pairs of vertices that are time-like separated on the original lattice implies that the 
intersection of any pair of red and blue domains enclosed by red and blue loops 
cannot be multiply connected. Similarly the absence of double triangles from the 
original simplicial configuration means that the (one-dimensional) intersection of a 
double line of one colour with a given domain inside a loop of the opposite colour 
must be either empty or simply connected. 

An important consequence of these constraints is that the number N22{t) of dual 
four-valent vertices is constrained both from below and above in terms of A'^i3(t) 
and A^3i(t). This does not happen in the matrix model where these numbers 
are completely independent. - It is possible that some of the regularity conditions 
discussed here can be relaxed without affecting the universal properties, but for the 
3d Lorentzian model this question has not yet been explored. 



Appendix 3 

The numerical simulations presented in this paper were performed for system sizes 
of 4k, 8k, 16k, 32k and 64k tetrahedra, and for total proper times T=16, 32 and 64. 
As usual, the standard unit was taken to be one sweep of the system, interpreted as 
iVa attempted moves. Since the acceptance of moves is a function of kQ, in order that 
all moves were performed approximately the same number of times, we had to tune 
the number of attempted moves for each of the three types of moves appropriately. 
This technique has been applied successfully before in three and four-dimensional 
simulations of Euclidean dynamical triangulations. We considered gravitational cou- 
plings in the range between fco = 2.0 and fco = 7.0. In this range the acceptance of 
move 1 is between 13.0% and 75%, the acceptance of the moves 2&3 between 3% 
and 11% and that of the moves 4&5 between 16% and 10%. For even smaller ko the 
acceptance of the moves 2&3 decreases rapidly and it becomes very difficult with 
the present set-up to change the geometry of the spatial intersections. A typical 
run corresponded to 10^ sweeps at a given value of ko. For all measured quantities 
we found autocorrelation times below 100 sweeps, which was also the typical time 
between successive measurements. 
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